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SECTION I 


INTRODUCTION 


This is a computer program which simulates air-to-ground gun 
effectiveness against a stationary target. The impact points of 
the individual rounds are correlated, and a Monte Carlo method is 

. required to determine the probability of hitting a rectangular target 
with one or more rounds in a single burst. The model assumes that 
gunnery is a stationary Markov process and that the aiming and 

t ballistic @ispersions are independent in the range (along the flight 
path of the aircraft) and deflection (normal to the flight path in 
the horizontal plane) coordinates. The guns are assumed to be fixed 
(2s opposed to guns turreted). The strafing aircraft flies at a 
constant airspeed and dive angle from a specifiicd slant range. 
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SECTION IT 


GENERAL DESCRIPTION : 


This 4s a Monte Carlo simulation program which determines the proba- 


bility of killing a ground target from an aircraft equipped with a : 

machine gun firing a single burst of N rounds. The individual aimpoints oe 
cannot be specified in advance; therefore, the program assumes that ; 4 
successive aimpoints are correlated (Reference 2). The aimpoints are i 


normally distributed about the center of the target. ‘he program further : 
assumes that ballistic dispersion is present. Thus, the 4*h round impacts 


not at its aimpoint but at some point nearby. 


The target assumed in the program is the rectanguiar projection of 
the real target on the plane normal to the line of flight of the attack- 
ing aircraft. Distances are measured with respect to a range-deflection : 
(R, D) coordinate system. The origin is located at the center of the i 
target, and the coordinate axes are parallel to the sides of the target : 
(Figwes 1 and 2). In real life the target will seldom or never be 
recta-gular in shape, but the projected target can always be approximated 
by a rectangle. The target length, 2, is measured in the range direction, 
and width, w, in the deflection direction. The slant range is denoted 
by s, c is the aircraft speed, and R is firing rate in rounds per minute. 
The slant range decreases for each successive round fired, thus increasing 
the apparent size of the t-rget and causing a corresponding inerease in 
hit probability. In addition to the above, the program inputs include the 
maximum number of rounds per pass (FN), the number of Monte Carlo itera- 
tions to be made (F), the burst length print increment (DN), the maximum 
standard deviation of the mean (E), the probability that the gun jams 


(PJAM), the number of gun systems (GUNS), and a time-to-rate table for 
a Gatling gun effect. 


To. determine aimpoint error, lev (Ry, Dy) be the aimpoint of the ith 
round, and let (ry, dy) be the point at which the round impacts. Further, 
let sq be the slant range at the instant the 1th round is fired. Then 
from Figure 2 it can be seen Ry = sqaq and Dy = 8484; expressing the angles 
in mils, Ry 0.001 syay; and Dy = 0.001 s4f4. The program uses the two 
angles a4 and 84 rather than Ry and Dy. These angles are computed by 
employing a random number generator. The angle o4 is @ random normal 
variable with mean zero and standard deviation og, and 64 is a random normal 
variable with a mean zero and standard deviation og. Bsllistic dispersion 
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4s cetermined similarly. Two angles y4 and 64 are defined so that (rq - Fy 

*9 001 8474) and (dy - Dy) = 0.001 sq4y. The angle y; 1s a normally 
Atecributed random variable with mean zero and standari deviation o,; 

8s, is a normally distributed random variable with mean zero and standard 


The range component of error is considered to be independent 


deviation Og: 
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Figure 1. 


Slant Range Geare try 


Figure 2. 
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of the deflection component. Also, for all i, j SH 


ECoya;) = alt - J\g2 
and 
E(B,8)) = pit Joa 


These two relationships characterize the programs correlation feature. 
The numbers a and b are inputs that determine the correlation between 
aimpoints. The correlation in range between consecutive aimpoints is 
given by a, and b -. Pos seta am deflection Bi genes in 


range between the ae ei the aimpoint is ali and the correlation 
in a tion is b . a=b#1, then Ry = and Dy = 94 for 
alli, j -%. This ce in the salvo model where every round has the 
sane aimpoint. If azb=0, every aimoint is independent of every 


other almpoint. The eae case will be when a and b lie between 0 and 1. 
It should be noted that if |1 - j| becomes larg, 


alt - 3! 45 anapli-Jl+o 


and the correlation Camps out. 


The correlation coefficient in range is a, and b is in deflection. 
These values can be inrut into the program if they are knovm. If the 
correlation coefficients are not known for the particular gun system 


being analyzed, the program will calculate them based cn the following 
equations: 
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where T is time in seconds between rounds. 


In the computation of the target kill probabilities, each iteration 
of the Monte Carlo process represents one pass at the tarmet. The 
attacking aircraft commences firing at the start of every pass and con- 
tinues until one of the following occurs: (1) a killing hit has been 
scored, (2) the gun jams, or (3) N rounds have been fired. When a gun 
jams, the number of guns in the system is reduced by one. ‘Therefore, 
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there is a corresponding reduction in rate of fire. If a Jam occurs, 
the firing loop is exited and a miss is recorded for calculation 
purposes. Zach round has a conditional kill probability (Py). This is 
the probability that a hit kills and it has the same value for every 
round in the burst if only one conditional kill value is input. A con- 
ditional kill value can be input for the beginning of the firing run and 
one at the end. If this option is elected, the program will do a linear 
interpolation between the beginning and end conditional kill values for 
each round based on tine into the burst. The program computes an ainpoint 
and a point of impact for each round and assesses damage by testing to 
see whether the round hit the target and if so whether the hit resulted 
in a kill. 


The final probability of kill is set equal to the number of 
successful passes divided by the total number of passes. 
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SECTION III 
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: 
I : a. MATHEMATICAL PROCEDURES 


The final solution is obtained as shown in the flow chari and 
the following mathematical procedures. Before starting the M-st 
iteration for the first data set dummy passes are made through the 
j random nurber generator. 
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For each Monte Carlo iteration there is given the standard devia- 
tion of the aim error (op, op) and the aimpoint of the (n - 1)st round. 
| There are also two Gaussian-distributed (mean = 0, standard deviation 
= 1) random numbers (a, y) selected, and the aimpoint of the nth round 
i (Rn, Dn) is determined by the following relations: 


for the first round (n = 1), then 


Ry, = oRA1 Cj) 
| dD * SnYy ; (2) 
ec a0 

Py = BR t ORC) nyt (3) 


2.% 
Dy * ODay + OplOO) yy (4) 


where a and b are the comelation coefficients in range and deflection 
and Ry, and Im are the respective range and deflection almpoints for the 
nth round. 
Equation (3) may be rewritten so that 
n 
4 -_ 
R, = On fav toy + (lea®)* Day (5) 
R 402 1 


and Equation (4) may be rewritten as 
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Then it can be seen from Equations (5) and (6) that 


E(R,) = E(D,) = 0 (1) 
E(R,°) = on 2 3 (8) 
B(D,*) = up” (9) 
B(R.R.) = op¢a Imm (10) 
E(D.D_) = op-b Imm (12) 


Therefore we have a process where each aimpoint is normally distributed 
with mean zero and standard deviation orn, op. We also have p process 
where the correlation coefficient ae the nth and the mt" round is 
alf-m| 4n the range direction and b/"-™] in deflection. When a and b 
equal zero or one the process degenerates into the independent or salvo 
cases, respectively. 


Next, the target size must be considered. The target dimension 
(2, w) 1s given with the firing rate (R), the aircraft speed (C), and 
the initial slant range (S). From these the half-target size (in mils) 
is determined for the instant at which the nth rourd is Mred using 
relations: 


Yn, 8 ————-—__ (12) 
S - (n-1)c(R) (1.688) 


500w 


5 (13) 
S = (n-1)c(R) (1.688) 


V2u, = 


The ainmpoint for the nth round is then checked to determine if it is 
within three standard deviations of the ballistic dispersion (fp, Bp) 
measired from the target center, or stated mathematically, 
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[Fol £2/22, + 38R (14) 
[Dy] =2/2un + 38D (15) 


If the aimpoint is not within three standard deviations of the ballistic 
dispersion from the target in either coordinate, it is assumed that the 
round missed the target and the aimpolint for the (n + 1)st round is then 
computed. If the aimpoint is within three standard deviations from the 
target in both coordinates, a Gaussian-distributed random numer (6) 1s 
determined and checked to ascertain whether the round falls within the 
limits of the target in the range coordinates. 


[Fn + BRI <1/2%y (16) 


If the round does not fall within the target limits a miss 1s assumed, 
and the aimpoint for the (n + 1)st round 1s computed. Hewever, if the 
nth round does fall within the target limits in the range coordinates 
another Gaussian-distributed random numer (E) 1s selected and the 
impact point of the round in the deflection coordinate is determined. 
The following check is then made to see if the round hit the target. 


[Dy + BpE] <1/2uy, (17) 


If the rouné does not hit the target, the almpoint for the (n + 1)st 
round 1s det-rndned. But, if the nth round does hit the target, a 
uniformly  ‘ributed (between 0 and 1) random number (PP) is selected 
and compar 1th the conditional kill probability to determine if the 
hit results in a target kill. If <PP, the weapon does not destroy 
the target and the aimpoint for the (n + 1)st round is determined. But 
4f Puy 2PP, the n&D weapon does kill the target. The count (Ny) of 
the e ‘arlo iterations for which the target is destroyed is 
increased ty one. Also, a counter (JJ,) for the next highest multiple 
(1) of the increment (AN) in the number of rounds for which the nroba- 
bility is to be determined is increased by one. 


M+ 1 +My , JI, +1 + 324 


When a round has resulted in a kill or when the maximm number of 
rounds (N) has been fired without killing the target, the entire precess 
is repeated until (F) Monte Carlo iterations have been completed. 
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Finally, the estimated probability of destro t 
N rounds is determined by: : ying the target with 


Nh 
Sr (18) 


The probability of destroying the target with J rounds is determined by: 
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SECTION IV 
PROGRAM UTILIZATION GUIDE 


The utilization guide for the air-to-surface gun simulation ccm 
puter program is contained in this section. The program varlables 
used are discussed with input formats, limits, and units specified for 
ean. 


Throughout the utilization discussion, variables which begin with 
the letters I, J, K, L, M, N are integer values and are right adjusted 
in their specified fields with no decimal punched unless otherwise 
specified. The alphanumeric formats are designated in the description 
colum of the program set-up procedure. All other variables are in 


decimal or real mode and may be punched anywhere in their columar fleld 
with a decinal point. 


This program has three different set-up procedures of which two will 
acconplish the same end result. The third set-up is a plot option only. 


Tne first set-up is a regular mm with any nunber of cases desired, 
end limited only by computer time. The second set-up shows the procedure 
for generating the input data when large parametric runs are needed. 

The third set-up describes the inputs for the plot only portion of the 
programs and should not be confused with the first two set-ups which also 
produce plots as desired. 


Figure 3 1s a flow chart overview of the gun simu.ation program 
with notes. 


The regular computer set-up is very sinple and straightforward 
but the generate-the-input set-up can be very ambiguous and fiustrating. 
When the second option can be used it should yield a 10 to 1 savings in 
set-up time. The analyst should become familiar with this option in 
order to better utilize his tine on large parametric computer runs and 
data requirement deadlines. It 1s hoped that the following flow chart 
will further simplify the generate-the-input set-up. 


TIMING 

The program run time is based on a number of variables in the 
program. It 1s not possible to figure the exact time required for 
each run. When the analyst becomes familiar with the program it wil). 


be easier to estimate the required mm time based on the following 
equation: 


RUN TIME = MCI ®N ® (NO. OF CASES) * .000159 
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where MCI is the number of Monte Carlo iterations 
N is numer of rounds fired on each burst. 
After calculating the estimated run time for a given set-up it is 
recommended that more time be added to the job card if the run time 
is not a limiting factor on the computer system. 


INPUT/OUTPUT 


Input and output are discussed in detail in Sections V, VI, and 
i VII. @ 
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OVERLAY 0,0 NOTE: Main Overlay Driver 


TTEST = 1 { NOTE: When ITEST=]1 it will be a 
. regular run 


A PERLA NOTE: Overlay 1,0 computes 
the results for the program 
NOTE: This decision determines 


whether the data has to be 
generated or not. 


NOTE: This system sub- 
routine prints out the 
Anput and generated data 


NOTE: The first data card 
is read which contains the 
program options. 


NOTE: At this point if the . 
set up is a regular mn it 
will begin computation 


Figure 3. Overview Flow Chart for dun Simulation Program 
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H NOTE: When ITEST=2 it 

: i will be a generate the 

ir ! input run 

i ae 

iy NOTE: Control has been RETURN TO 

: returned to the driver OVERLAY 0,0 

t (Overlay 0,0) 

i : 

t | NOTE: Overlay 2,0 will 
! generate the input data 
bo 

t | 

F NOTE: First read tape 5, 
t then tape 4 
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NOTE: ‘The decision here is 
whether to execute the 
generated input or not. 
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SECTION V 
REGULAR COMPUTER RUN 


SAMPLE PROBLEM 


The effectiveness of a strafing tactic for a 150-knot aircraft with 
a firing rate of 4100 rounds per minute is to be analyzed. The pilot 
commences firing at a slant range of 1000 feet and fires a single burst 
of 120 rounds at an 11.74 x 11.74-foot target. The standard deviations 
of aiming error are 4.24 mils in range and deflection, and the standard 
deviations of ballistic dispersion are 1.39 mils in range and deflection. 
The probability of kill given a hit for the first case is a three-round 
mix with the ammnition belt having four rounds with a PHK of 0.012, 
0.016, and two rounds with 0.020, 0.024 and one round with 0.006, 0.010. 
There are two probability of kill values in each case. ‘The first value 
is the probability of kill given a hit for the beginning of the burst, 
ana the second value is for the end of the burst. The computer does a 


linear interpolation between these two values based on time into the burst. 


Table 1 contains the description for setting up any regular computer 
run. Table 2 1s a sample set-up of the case described above. Table 3 
is an output listing of the input data. Table 4 contains the final out- 
put probability of kill for only two of the cases. The first case was 
at 1000 feet slant range and 150 mots. The second case shown in Table 4 
was for 1000 feet slant range and 350 knots air speed. The second case 
has deen included to show the analyst what happens when the conditions 
become unrealistic. This feature saves paper and computer time. Ficure 
4 contains the probability of kill for the first six slant ranges. Card 
2 in Table 1 indicates the number of cases plotted on one graph. 


The required time to run any given set-up is described in Section IV. 
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TABLE 1. DESCRIPTION OF THE SET-UP FOR A REGULAR COMPUTER RUN 


Descrintisn 


Variabie 


DESCRIPTION OF TI SET-UP FOR A REGULAR COMPUTTR RU? 


1 IPLOT 0,1.2 0 = tl plots 
1 = Plot option has been turned mm 
2 = Stand alone plot (see PLOT OILY SET-U) 
IOPT 0,1 0 = Regular computer run, leave IGCO Liark 
1 = Data will be generated (see GITRATS 
TIE INPUT SET-UP) 
1Go 0,1 CG = Do not execute the generated incut 
data 
1 = Execute the venerated input data 
(see GENERATE THE INPUT DATA SET-UP) 
‘ITS 1-10 Number of burst lengths to be plotted 
IPLOT(1) Time which burst leneth should be nlottec decimal 
IPLOT(2) 
69-75 1PLOT(10) Tenth burst length time decimal 


NOTE: Card 1 is aeressary in the set-up even when all the options are 
7070. 


1-5 SLAIT 1-5 If IPLOT = 1, you must input the number 
of slant ranges to be plotted on one 


raph. If no plots ar: wanted, this card 
must be omitted from the set-up. 


3 1-60 TITiF. tle or general information (alpha- 
numric data) 
61-70 DIVE £90 Dive mele - for identification miv NES 
4 1-2 D1) 1 Address = 1 . 
3-10 A 20, $1 Correlation coefficient in renee beteen 


consecutive almpoints. 


NOTE: A value for the correlation coefficient is computed 


in the procran for each round. This parameter should 
be set to zero {f JBI = 1 on cari 2h. 


iad. weak bleed Sas hii at 


TABLE 1. 


DESCRIPTION OF THE SET-UP FOR A RHOGULAR COMPUTER RUN (CONTINUED) 


fg me 


6 A 8 REP AH 
eel 


Address * 7 


j : Slant range at the beginning of 
bag firing run. 


Address = 8 


ea TRE, 7 ee ET 


£ ‘ 3-10 Firing rate of gun ir: rounds per Pe 
i : minute. -3 
i ' le 1-2 Address = 3 1] 
t 

F 3-19 Aircraft speed mots q 
t 

b 13 1-2 Address = 10 { 
' * 
é 3-10 Number of rounds Mred on a single decinal j 
} pass per gun. 3 
¢ 4 
: 1 1-2 0611) li Address = 11 

3-10 NTYPE 1,2;3 Nunber of tynes of mixed belts. Set 


equal to 1 if no mixed belts. 


\ 


yaad li tus nl A al ii a 


15 1-5 NUMR(1) Blank-N Number of consecutive rounds using 


this conditional kill probab‘lity. 
If left tlank the prepram sets 
NUMR(T) equal to N. 


Starting value for the conditional 
kil] probability - 4f CP! ts left 
blank, CP1(I) is used for ali rounds. 


‘ 


Tul Poor oa ua bilan. ute Basan Ste ans ll 


Fnd value for the conditional kill 
probatility. If an end value is used, 
the program does a linear interpolation 
between CPl and CPN based on tire. 


POTS TN RRR PT ne AE ST CHT eT ATaR STS 
Nn 
J 
= 
WwW 
rary 
~ 
2 
a 


NOTE: Repeat card 15 NTYPE times. 


1-2 D(12) 12 Address = 12 

3-10 i Target length feet 
decimal 

1-2 D(13) 13 Address ® 13 

3-10 W Target width feet 


t 


aan! i cL a 
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TABLE 1. DESCRIPTION OF THE SET-UP FOR A REGULAR COMPUTER RUN (CONTINUED) 


H Descriptic:. 
E : 
t ar) Address = 14 ae 
.7 
E s 
4 . Fj 1 
t Maximum number of Mente Carlo iter~ i 
f at ions. : 


4 
Hy 


Address = 15 


Number of dummy passes ihrouph random 


decinal 
numer generator, 


Address = 16 


Increment in durst lensth. 


Tris controls the number of lines that will be printe:i; i.e., 1 


will cause the nrinter to write out a data line for each round 
in the burst. 


al 1+2 ™17) 7 Address = 17 
3-10 c Desired maximum value of tie standard 
deviation of the mean, 
22 1-2 g -D(18) 13 Address = 18 
3-10 PIAM ~obability of the mun Jameng. 
33 1-2 ™ 19) 19 r_wess = 13 


3-19 GUS Nusber of mun svstems .o te analyzed. 


NOTE: The orozran computes a final probability of kil] based an the 
total number of gin systems. 
24 1-2 (29) 20 Address = 20 


3-16 JIM G41 


= indicates that vou have input som 
correlation value in cards 4 and 5 
other than zero and omit cards 25, 26 
and 27, 

1 = indicates you have a zero in cards 
4 and 5 and plan to input a time-to~rate 
table for a Satiing mm Ly completing 
cards 25, 26 and 27, 


decimal 


tt an faa AR A ae at te a al 


1-5 NOT SG Number of pairs of entries for the time- 
to-rate table for a Gatlin® mun. 
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TABLE 1. 


26 


27 


28 


Number of rounds fired at TIME (1) 
Number of rounds fired at TIME (2) 


. 


Number of rounds fired at TIME (NOT) 


NOTE: If JIM = 1 you can have 1S cards £3, typical for card 27 also. 


1-5 TIME(1) Time to fire RD(1) rounds 
6-10 TIME(2) Time to fire RD(2) rounds 


Time to fire RD(NOT) rounds 
BLANK CARD AT END OF EACH DATA SEY 


3 cards | TIME(NOT) 


NOTE: For multiple cases, include the title (card 4) and any 
parameters that may change from card 4 to 28. Cases 
ars untinited. If the time to rate tazle (cards 25, 
2c, 27) is changed, tue conditional kill probability 
(cards :4, 15) must also ve repeated, 
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DESCRIPTION OF THE SET-UP FOR A REGULAR COMPUTER RUN (CONCLUDED) 
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SECTION VI 
GENERATE THE INPUT RUN 


SAMPLE PROBLEM 


The sample for Generate the Input Run uses the same data as the 
Regular Computer Run in order to demonstrate the flexibility, time 
savings and less chance for errors that can be realized through the use 


of this technique (see sample problem in Section V for details of problem). 


There are five parameters that change more frequently than the other 
inputs. They are listed in Figure 5 with their hierarchy. Table 5 
describes how to set the job up. Table 6 is a sample coding for the 
problem. Table 7 is a REPRNT listing (a computer system routine) of the 
input data. Table 8 contains a list of the data generated “from the input 
from Table 7. The final output Misting is given in Pale 9 with most of 
the input values included. After the columar t4tles each round is 
analyzed for the entire burst. This output +4@ the same as the Regular 
Computer Run. Also when the slant range’ oes negative the computation 
stops and prints out a message thatthe aircraft flew into the target. 
Figure 6 is the filmplot output~for the first six slant range values with 
their respective probability of kill values plotted. 


a 
The time requtred for any given computer run can be calculated by the 
equation described in Section IV. _ 
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AIM 
ERROR (S) 


1 
DO 6 ,I=1,LAIM 
TAIM >1, <3 


oo abate stk let 


BALLISTIC 
DISPERSION (S) 
IBALE >1, <6 


2 
DO 6 ,J=1, IBALE 


AIRCRAFT 
VELOCITY IES) 
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DO 6 ,K=1,IVEL 


IVEL > 1, < 12 

TARGET (S) ; 

LENGHH § WIDIH : 

ITGTL ' 

ITGIW j 

SLANT 5 

RANGES DO 6 ,M=1,NSLT(6) 

NSLANT >1, < 6 

Figure 5. Flow Chart for Generate the Input Run 
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TABLE 5. 


DESCRIPTION OF THE SET-UP FOR A GENERATE THE INPUT RUN 


fescription 


Yo 


Thr. 


0 © No plots. 
1 = Plot option has been turmmei on, 
2 = Stand alone plot. 


3 = Regular computer run. 
] 3 Lata wlll he mnerate). 


pat 
ae 


9 = ho net execute the generated input 
juta. This ts used fer checking 


oyeteus before excarineg It. 
Tf ocaeonte the generates ioput vata. 
vanes 1-20 terteor of burst leneths te be plotted. 
Type TEL) Fine ahton burct leneth choald be plotted, 
rem?) 
rey 19} OUR ME), 
vA 


secs 
decimal 


ney 3 Addners 31. 


3-12 4 25,8) Sormelatior coerMetent te rancr 
between consecutive alrpotnts. 


‘it: Tf a Tw for the correlation coefficient is camputed in th program for 
eagh coun}, his raramter should be cet to vem. Af JT = i oon cari 12. 


1-2 a) 2 Address = >. 


215 RB 2, $1 Cormelation coefMelent In aioMect ton 


h ion consecutive aimpolntc. 


ate: Tf a value for the correlation coefficient is computed in the 


rrogram for 
cach, round, “htc parareter should te set to gerne if 704 
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DESCRIPTION OF THE SET-UP FOR A GENERATE THE INPUT RUN (CONTINUED) 


Description 


Address = 8. 
Firlig rate of gun in rounds per minute. 


Address = 10. 


Number of rounds fired cn a sinzle pass decimal 
per gun. 


Address = 14, 
Maximum number of Monte Carlo iterations. 
Address = 15. 


Number of dummy passes through random decimal 
number generator. 


Address ® 16. 


Increrent in burst length. 


NOTE: This controls the numer of limes that will be printed, 1.e., 1 will 
cause the printer to write out a data line for each round in the 


3-10 
1-2 


3-10 


burst. 
D(17) 
E 


D(18) 
PJAM 
D(19) 


GIS 


17 


18 


19 


Address = 17. 


Lesired maxirum value of thre standard 
deviation of the mean. 


Address * 18. 
Probability of the gun janming. 
Address = 19. 


Numer of gun systems to be analyzed. 


NOTE: The program computes a Sinal probability of kill based on the total 


nunber of gun systers. 
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e eo TABLE 5. DESCRIPTION OF THE SET-UP FOR A GENERATE THE INPUT RUN (CONTINUED) 
r 
is a 
E Card Colums | Varlable Limits Deccript lon Units 
c: 
b 12 1-2 D(20) 20 Address = 20 
} [ 3-10 um 0,) 0 indicates that you have input come 
f f correlation value in cards 2 and 3 
E { other than zero and omit cards 23, 14, 
} t and 15, 
; } 1 indicates you have @ zero in camis 
t y 2 and 3 for tie correlation values und 
7 E plan to input a time-toerate table for 
7 | a Gatling gun by completing cards 13, 
; t 14, ord 1b. 
; } 
: t 13 1+, Ww S¥. Number of pairs of entries for the 
t | time-to-rate table for 1 Gatling sun. 
i 1 1-8 HC) hurber of rounds Med at TRE 0) 
' i 9-16 R23} Surber of rounds Mined at TIME (ly 
} i ' 
{ i 
i ? cords | RO(NCT) Number of rounds Mined at TIME (HUT) 
: } WOM: You can have from 1 to 3 cards for a myimin cf BD entries. 
z : 
t i 1g 1-5 TIKES 1) Tim to fine KOC1) round: 
t i Geld | TIME?) Mee te thre BCE) round: 
f. i i 4 : 
3 cards | TIME(MUT) Time to Sire RO(NUT) rewsac 
b t 
; : 6 Blank card 
i ry End of record cam 
f- : 
& : 
4 i 13 1-2 TAIM Lp4 Number of acts cf alming ermro 
t i 6~10 STGIAL Beginning value for the frat otundard mile a 
E { (ALPHA) /L deviation of the aimerror in range 
§ : 21-35 GI) fA valu for standard deviat¢on alm mis o 
i : error Tn range 
| NOE: Tf no end value de used, the beginning value will be used for the - 
: : ertie: burst leqmh. Tf an ent value ie use), the provran. will 
doa ‘Miner tntverveiation tetacen Gs veginning aviend valu for 
Uae Abn nor Laced on the aloe’sy canes of the afreraty, 
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TABLE “. DESCRIPTION OF THE SET-UF FOR A GENERATE THE INPUT RUN (CONTINUED) 


ae am 
16-20 SIGH Teeinnine value for the first standard mils o 
(BETA) (1) deviation ef the air error in deflece 
tion. 
21-28 3GD1L(1) Tnd value for standard G@viation aim mils v 
error in deflection. (See SGRL(1) Note) 
26-H STIGMA Beeinning value for the second standard 
* ALPHA) (2) deviation, of the aimcrror in rane 
31-35 S521(2) Dui value for the seconds standard 
deviation ain error in renga. (See 
go. S'PI.61) Mote) 
36-40 STW Repinninre value for the second standard 
(BETA) (2) deviation alm error in deflection. 
G1-hs $4D1(2) nd value for the second stendard 
deviation aim arror in deflection. 
(3ae SGR1(1) Note) 
46-50 SIGUA Bepinnine value for the third standard 
deviation aim ervor in rence 
$1<£$ §GR1(3) Ind value for the third standaru 
deviation aim error in rence, 
(See SGHi(1) Note) 
56-60 STIGMA Beginning value for the third standard 
(BETA) (3) deviation aim error in deflection. 
61-65 $4D1(3) End value for the third standant 
daviation aim orror in deflection, 
(See SGRL(1) Note) 
19 1-2 TBALT 1-6 Number of sets of ballistic errors 
€~-10 SALE. 2) First ballistle error standard deviation | mils o 
rane 
11-18 BALE(2) Tirat ballistic error standard deviation | mile 9 
dn deflection, Repest the above stop 
and this step for a madmun of 6 sets. 
90 1-2 IVEL 1-12 


Number of aircraft velocities following 
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TABLE 5, DESCRIPTION OF THE SET-'IP FOR A GENERATE THE INPUT RUN (CONTINUED) 


22 


% description 
6-19 VEL(L) First aircraft velocity i 
VEL(12) Twelfth aireraft velocity | ants 


TTTL Number of different targe’. leagins 


following ; ae 
TOTL(1) Pirst target length in range 7 tye 
Twelfth target lengh in range yoo ft feet 


Number of target uldths following 


Mret tanet width in deficcti in feet 
fecina. 

615 | ITGIw(42) Twelfth target width din deflet Lon feet 
a~2 TKI Number of conditional Kila caras tn | deci? 


this set - must equal NGLT(1), Card 
eu 


NOE: This value is used on the Mrst card of each set of conditional 
Will ratio tables, t.e., whan card 23 1s repeated second through 
the sixth time thic parameter io omitted. See sample sect up. 


6-19 NTYFE‘1) ae. l,e, or 2 Jiffer-nt arm types analyzed decisnal 
in burst 
1-15 aK) Nurber of consecutive rounds using the decimal 


first conditional kill ratio 


WT: Tf MUMR(L) = 6 the progree sets ‘UMR(1) equal to if MVYPE(1) #1. 
This card typ: is the same as cand 24 in the REGULAR COMPUTER run 
BCL. 


14-20 Cha ees Value of the conditional kill ratio at 
the beginning of the burst 
21-25 CPIL1) 


Value of the conditional kiil ratio at 
the end of the burat. The program does a 
linear interpolation between the begin- 
ning ami end values. If no end condi- 
tional kill value is irmut, the hegin- 
ning value will be used MUMR(1L) times. 


EE A AK A ET A TR RT ANS, FT NETO Oe TEN ET TE LE I TT I A TE TT 


| 
! 
| 


Sng dasa eae ee 
dl 
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= 


NUMR(2) 
31-35 cP1(2) 
BNO CPN (2) 
41-05 NUMR (3) 
46-50 CP1(3) 
51-55 CPN( 3) 


Number of consecutive rounds using the decimal 
second conditional kill ratio. This 

opticn will be used in NTYPE(1) is 2 

or 3. 


Beginning value for second conditional 
kill ratio in tne mixed round ammo belt. 


End value for the second conditional ki11 
ratio in the mixed round anmo belt. When 
the end value ia input, the program does 
a linear interpolation between the begin- 
ning and end values. If no end condi- 
tional kill value {s input the begiming 
value will be used NUMR(2) times. 


Number of consecutive rounds using the decimal - 
thira conditi.:rsl kill ratio. This ‘ 
option will be used if NTYPE(1) is equal 

to 2. 


Beginning value for the third condi- 
tional kill ratio in the mixed round 
ammo belt. 


End value for the third conditional kill : 
ratio in the mixed round amm belt. When “ 
the end value is input, the program does ad 
a linear interpolation between the begin . 
ning and end values. 


NOTE: Repeat Card 23 NSLT(1) times for each additional set of condittorial 
kill tables (Sets range from 2 to 2ulu). A set 4s = NSLT(1) or 
IKILL. No. of sets required = ITAIM*IBALE®IVELFITGTL. 


24 1-5 TSLAG 1 Nurber of unique slant range tables; see 
sample set up. ‘ 
6-10 NSLT(1) 2-6 Number of slant ranges on this card. 
NOTE: NSLT(1) should equal IKILL. 
11-20 SLRNG(1) Mirst slant range : feet 
61-70 SLING (IKILL) slant range : feet 
(IKGLL) : 
PA 1-60 TITLE(1) Hollerith information 
61-70 Crit this parameter on GENERATE~THE- 


INPUT data set~up. 
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SECTION VII 


STAND ALONE PLOT 


There may be times when the probabtlity of kill values are known 
but no plots were made. This section becomes a back-up plot option when 
this occurs. The Pe's can be keypunched and plotted with a minimum 
computer and turn around times. Without this option ths analyst would 
have to resubmit the computer run with the plot option tumed on or plot 
the Px's by hand. Table 10 describes the set-up and Table 11 is a 
sample set-up. Figure 7 1s the filngiot output for the stand alone plot. 
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APPENDIX A 


FORTRAN VARIABLE LIST 


This appendix contains a list of input, output, and intermediate 
FORTRAN variables used in the mathematical computation overlay 1,0. 
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FORTRAN 


Variable 


~ A,DCL) 
AT 


ALPHA 


B,D(2) PDC) 
PBB , TEMP 
BBC 


BETAD,D(6) 
BETAR,D(5) 


BURSTL 
C,D¢9) 
CPK(T) 


COUNT ,12 


CPKI “T) 
CPK2(I) 
CPK3(T) 


(71(T) 


CPN (I) 
mr 


DELTA 


Description 


Correlation, coefficient in range - input variable. 


4n intermediate variable used to check to see if the time- 
to-rate table is correct. 


A xandam normal number used to determine the aim point 
in renge. 


Intermediate aim error value in range. 

Intermediate aim error value in deflection. 
Correlation coefficient jn deflection ~ input variable. 
Intemnediate aim error value in range. 

Intermediate aim error value in deflection, 


Standard deviation of ballistic error in deflection - 
input variable. 


Standard deviation of ballistic error in range - input 
variable. 


Burst length input variable for plot oly opticn. 


Input and output variable for the aircraft speed in Imots. 


Conditional kill probability, or probability that a hit 
kills - output variable. 


The number of times the program goes through the Monte 
Carlo loop. 


lst conditional kill value using a mixed anmw belt. 
2nd conditional kill value using a mixed ammo belt. 
3rd conditional kill value using a mixed ammo belt. 


Starting conditimal kill value at the beginning of firing 
run. 


End oonditional kill valve at the end of firing run. 
Aim point of the first round in deflection. 


Randem normal number used to check the ballistic error 
in range. 
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FORTRAN 

—Vartable OL eseription 4 

DIVE Dive angle (for information only), input, and output . 4 

variable. | 

DN ,D(16) Output print increment in burst length - input variable. E 

ne : An intermediate value used to calculate time into burst 4 

for each round. 3 

or Time increment between rounds in a burst. ’ 

: DUMMY Acts as a return variable to mix the nandom number generator 4 

-- not used in computation. j 

E,D(17) Maximum allowable error for standard deviation of the mean- F 

input variable. 3 

EPS Random normal number used to check the ballistic error F 

; in deflection. | 

F,DC14) Maximum number of Monte Carlo iteratians ~- input and output 4 

variable. | 

FF(I) Output probability of kill for each round (AN). 4 

FI (12) Target length in feet-input and output variable. i 

FIN2 ,FLR2(I) Half of target length in mils. 4 

“S 

FLR2(1)=TEMF Half of target length in mils. 4 

FN ,D(10) Attempted number of rounds fired om a single pass - input 4 

variable. q 

F 

GAMMA Random normal number used to determine the aimpoint in 5 

deflaction. | 

GUNN Intermediate variable used to increment gm system. F 

GUNNS ,D(19) Number of gun systems considered - input variable. 

IAIM Number of sets of aiming errors - input variable. 4 

4 

Ii Initial value of the index for the Monte Carlo loop, f 

12,0 Test value of the index for the Monte Carlo loop (2209). F 

; Each iteration represents one pass at the target. : 
| ION , IN Intermediate variable for burst length increment. 
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LST, INC 


NCASES 
NEMP ,D(15) 


NXILLS 
NOT 


NUMR(T) 


P,X11) 
PD(I) 
PJAM,D(18) 
PK 
PKPLOT(T) 


Description 


The initial and index values for the DO loop that deals 
with the 3 belt mix, or 3 conditional kill values for 
1 pass. 


fn integer value calculated for each round in the burst. 
Integer value of rounds per pass per gun - output variable. 
Integer value of the address on the input cards. . 


Integer value used to determine the number of lines to 
be printed on a page. 


Integer used in an intermediate calculation of the condi- 
tional kill probabilities of a mixed belt set up. 


Intermediate integer used with a mixed belt conditional 
kill probability set up. 


Attempted number of rounds fired on a single pass. 
Number of cases to be plotted. 


Number of empty passes through random number generator - 
imput and output varlable. 


Number of kills or successful passes. 

Number of paire of entries in the time-to-rate table. 
The number of slant ranges to be plotted on one graph. 
Number of burst lengths to be plotted. 

The mamber of sypes of mixed belts. 


The number of oconsecutive rounds using this conditional 
kill prebability. 


Input conditional kill prebability 

Correlation coefficient array in deflection. 
Probability that the gun jams - input variable. 
Probability of target kill - output variable. 
fn intermediate value used to plot the PK's. 
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RC 
RD(T) ,RPS 
RNSTL 


RPS 


RPS1(J), IR 


S,D(7) 
SDPK 


SGD1 


SGRL 


SIGD,DC4) 


SIGDUT) 


STGR,D(3) 


SIGRI(T) 


"The ena paremeter on the input data card if it has one. 


Description 


The parameter after the address on each input data card. 


A random number used to test against the conditional kill 
value to see if the target was killed. 


Correlation coefficient array in range. 


Title for the plots for the plot only optic in the 
program. 


The updated value for the index of the Monte Carlo loop. 


Steady state firing rate per gun in rounds per min - input 
and output variable. 


_ 
' 
3 

4 
4 
i 

7 | 
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} 


Aim point of the first round in range. 
Rounds per second for stand alone plot. 


Total rounds fired - used on plot only option - input and 
output variable. 


Rounds fired in one second - input and output variable 
for plot option. 


The round number array for plotting. 
Slant renge - input and output 


Output variable for the standard deviation of the mean - 
used to test accuracy of probability of kill. 


End input value fer the standard devi.ation ain error in 
deflection. 


End input value for the standard deviation aim error in 
range. 


Standard deviation of aim error in deflection ~ 7 
Anput and output variable. 


fn intermediate array variable for the standard deviation 
aim error in deflection. 


Standard deviation of aim error: in range - input arid output 


An intermediate array variable for the standard deviation 
aim error in range. 
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FORTRAN 


SLR(TI) 


T(1) 


TCODE 


TD(T) 
TD 


TEMF 


TIME(T) 


VR 


w,D(13) 
WCOLE 


WN2 ,WR2(T) 
x 
xX 


Variable 


Description 


Output array variable for slant range at the time each 
round is fired. 


Intermediate time array variable. 


Target code for the plot option oly - input and output 
variable, 


Delta time array between rounds. 


TDCI) squared - this variable is used to find the range 
and deflection correlation coefficients. 


Intermediate variable which converts 1/2 the target 
length to mils. 


Intermediate variable which converts 1/2 the target 
width to mils. 


fn intermediate variable used in the look-up table. 
Time value array for the time-toerate table. 
Description of the output in 60 characters or less. 

An array that contains the burst lengths to be plotted. 


The variable where knets are converted to feet per 
second. 


An intermediate variable that determines the distance 
plane traveled between rounds fired. 


Target width in feet - input and output variable. 


Weapon code used as an input - output variable in the 
stand alone plot. 


Half of target width in mils. 
Probability that the gum jase. 


Random number used io test against X to see if the gm 
janmed, 
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OVERLAY 1,0, 


PROGRAM RUN 


DIMENSION D(20), TITLE(6), TPLOT(10), PKPLOT(10), 
RPS1(10) 


COMMON JJ(1000), FF(1000), SLR(1000), FLR2(1000), 
WR2(1000), CPK(1000), TIME(30), RD(30), TD(1000), 
PR(1000), PD(1000), T(1000), SIGR1(1600), 

SIGD1(1000) , CPK1(1000), CPK2(1000), CPK3(1000), 
CP1(3), CPN(3), NUMR(3) 


COMMON/ENDDT/IEOF, ITEST, IGO, JIM, IPLOT, IOPT, 
| NTS, TPLOT 


EQUIVALENCE(D(1),A), (D(2),B), (D(3),SIGR), 
(D(4) ,SIGD), (D(5),BETAR), (D(6),BETAD), 
(D(7),S), (D(8),R), (D(9),C), (D(10) , FN) 
(D(11),P), (D(12),FL), (D(13),W), (D(14),F), 
(D(16},DN), (D(17),E), (D(18),PJAM), (D(19) ,GUNS) 


EXTERNAL TABL1V 


DATA TARGU/SHTY = / 


CALL PLTIDV 


Figure 8. Flow Chart of Overlay 1,0 
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ASSIGN 400 
TO IEOF 


SSeS ALS ie ea cen es a pe er en rr re SE ee ee ee 


OVERLAY 1,0] 


PROGRAM RUN 


DIMENSION D(20), TITLE(6), TPLOT(10), PKPLOT(10), 
RPS1(10) 


COMMON JJ(1000), FF(1000}, SLR(1000), FLR2(1000), 
WR2(1000), CPK(1000), TIME(30), RD(30), TD(1000), 
PR(1000), PD(1000), T(1000), SIGRI(1000), 
SIGD1(1000), CPK1(1000), CPK2(1000), CPK3(1000), 
CP1(3), CPN(3), NUMR(3) 


COMMON/ENDDT/IEOF, ITEST, IGO, JIM, IPLOT, IOPT, 
NTS, TPLOT 


EQUIVALENCE(D(1),A), (D(2),B), (D(3),SIGR), 
(D(4),SIGD), (D(5),BETAR), (0(6),BETAD) , 
(D(7),S), (D(8),R), (D(9),C), (D(30), FN) 
(D(11),P), (D(12),FL), (D(13),W), (D(14),F), 
(D(16),DN), (D(17),£), (D(18),PJAM), (D(19) ,GUNS) 


EXTERNAL TABL1V 


DATA TARGU/SHT4 = / 


CALL PLTIDV 


GO TO — 
(600 , 8} ITES 
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ASSIGN 400 
TO IFOF 


ECT TO RT TT agree Tre 


THIS OPTION 
DOES A STAND 
ALONE PLOT 

WHICH STARTS 
AT 1000 


READ IPLOT, 
IOPT, 1GO, NTS 
(TPLOT(I), I= 


CALL EXIT 


IF 
TOP?.FQ.1 & 
TPLOT.GT.1 


WRITE "INPUT 
CANNOT 8E 

GENERATED UNLESS 
PLOT = 1 or QO" 


CALL EXIT 


; Ficure 8. Flow Chart of Overlay 1,0 (Continued) 
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Flow Chart of Overlay 1,0 (Continued) 
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Figure 8. Flow Chart of Overlay 1,0 (Continued) 
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Figure 8. Flow Chart of Overlay 1,0 (Continued) 


72 


4 
4 


SMa sua La ta = le 


i 


Wie 
alin inal 


Ogee O ee ee he ge ee Ae ee eee ee gee ee pee, eee gy a ee 


<> (res) CPK1(1)=CP1(J) 
KS 


> (v3) CPK3(I)=CP1(J) 


Figure 8. Flow Chart of Overlay 1,0 (Continued) 
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JSTR = 1 
IST = 2 
C) DO IV=3, NOT 


(11s) JSTR=.JSTP 


IST=IST+1 


BBB=CP1(J )+(CPN(J) - 
CPL(J)) *(TIMEC IST -1)/ 
TIME (NOT) ) 


(ves BBB=TEMP 


Figure 8. Flow Chart of Overlay 1,i) (Con! tnued) 
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Figure 8. Flow Chart of Overlay 1,0 (Continued) 
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Figum 8. Flow Chart of Overlay 1,0 (Continued) 
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TIME (J-1))/(RD(J) 
-RD(J-1)) 


WRITE TIME 
VS ROUND 

NUMBER TABLY 
IS WRONG 


IF 


pelts lb sna 


seb 


DQe TIME (J) -DT#RD(J) 


Figure 8. Flow Chart of Overlay 1,0 (Continued) 
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Figure 8. Flow Quart of Overlay 1,0 (Continued) 
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Figure &. Flow Chart of Overlay 1,9 (Continued) 
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Figure 9. Flow Chart of Overlay 1,0 (Continued) 
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CALL RANG 
(ALPHA , GAMMA) 


RC=A*RCHSIGRI (I 
*SQRT(1.0-A*A) 
*ALPHA 


DC=B*DC+SIGD1(J) 
*SQRT(1.0-B*B) 
*GAMMA 


On 
@ 
CD 


PK=FLOAT(NKILLS) /COUNT 


SIMs AND 
SLR (JJJ)<0.0 


Figure 8. Flow Chart of Overlay 1,0 (Continued) 
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Figure 8. Flow Chart of Overlay 1,0 (Continued) 
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WRITE (6,187) 
SIGR,SGR1, 

BETAR,SIGD, 
SGD1, BETAD 


| 
; 
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IF 
JIM#0. AND 
SGR1<O 


WRITE (6,191) 
SIGR,BETAR, 
SIGD ,BETAD 


WRITE (6,200) 
¢ 


WRITE (6,205) c 
R rn 
WRITE (6,210) 
Ss 
KN=FN P 


WRITE (6,215) 
KN 


Flow Chart of Overlay 1,0 (Continued) 
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Figure 8. Flow Chart of Overlay 1,0 (Continued) 
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Figure 8. Flow Chart of Overlay 1,0 (Concluded) 
95 
Sao OTR mn ne carne a en OES SONI APM TNE Be teat SARA age aE Oe eR Nn pe Ty = 


i 
calla! 


= 


APPENDIX C 
P0655 PROGRAM LISTING 
‘ Appendix C contains a FORTRAN program listing of the air-to-ground 


gun simulation program complete with three overlays and three subroutines. 
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